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Abstract 

We consider two-level detectors coupled to a scalar field and moving on arbitrary trajec- 
tories in Minkowski space-time. We first derive a generic expression for the response function 
using a (novel) regularization procedure based on the Feynmann prescription that is explic- 
itly causal, and we compare it to other expressions used in the literature. We then use this 
expression to study analytically and numerically the time dependence of the response func- 
tion in various non-stationarity situations. We show that, generically, the response function 
decreases like a power in the detector's level spacing, E, for high E. It is only for stationary 
world-lines that the response function decays faster than any power-law, in keeping with the 
known exponential behavior for some stationary cases. Under some conditions the (time de- 
pendent) response function for a non-stationary world-line is well approximated by the value 
of the response function for a stationary world-line having the same instantaneous accelera- 
tion, torsion, and hyper-torsion. While we cannot offer general conditions for this to apply, 
we discuss special cases; in particular, the low energy limit for linear space trajectories. 



1 Introduction 



The physics of accelerated quantum systems in Minkowski space-times share some interesting 
features with quantum gravity phenomena that are much less tractable; e.g. with the Hawking 
evaporation of black holes[lJ and quantum aspects of expanding universes and inflation[2] . Beside 
being of great interest on their own right, such Minkowskian phenomena may shed light on their 
full hedged curved space time cousins. 

A moving two-level detector (so called a Unruh or de Witt detector) that couples to the 
modes of a field of interest is a simple and very useful paradigm for studying quantum effects 
of accelerated systems. The most renowned of the phenomena studied with such a device is the 
Unruh effect [3], whereby a linear, uniformly accelerated detector in flat space-time perceives the 
vacuum as a thermal bath of the modes to which it couples. (Whenever we refer in the paper to 
space characteristics of a world-line, such as its being linear, circular, etc, we mean it to be so 
in some Lorentz frame.) One starts by defining a vacuum-the ground state of the various fields 
of interest-based on a set of creation-annihilation operators associated with inertial observers. 
This means that modes of any free field are eigenmodes with respect to the proper-time of such 
inertial observers. 

The vacuum effects associated with the detector are encapsuled in the response function, 
which is the object we study in this paper. The correct way to regularize the divergent expression 
for the two-point Wightman function, which appears in the definition of the response function, 
is still a moot issue. For example, following work by Takagi[4J, Schlicht[5j proposed to start with 
a detector of finite spatial width in the proper frame. In this case, the Wightman function is 
regular, and the point detector response function is gotten by taking the width to zero. This gives 
well-defined expressions, which have been generalized by others [Sl[7j. We discuss the merits of 
this approach, and propose an alternative that is, in our view, superior in some regards. 

The response function is an explicit function of r, the proper-time along the world-line, 
and of E; it is also a functional of the world-line. Since the response function is Poincare 
invariant, it is reduced, in 3 + 1 dimensions, to a functional of the three world-line invariants: 
the acceleration a(r), the torsion T(r), and the hyper-torsion H{t). Beside the trivial, inertial 
case, only two world-lines have so far lent themselves to analytic evaluation of the response 
function: the original, Unruh case (constant a, T = H = 0) , and planar cusp motion (constant 
a = T, H = 0), which is also the infinite-Lorentz-factor limit of uniform, circular motion 
response function, as we show below. The high-7, circular Unruh effect has been discussed 
extensively (e.g., by Bell and Leinaas[8j, and by Unruh[9|), who obtained analytic results for the 
infinite Lorentz factor case, in the hope that it can be detected in electrons moving in storage 
rings. Letaw and Pfautsch|10| considered, numerically, circular motion with finite velocities. 
Experimental manifestations of the linear Unruh effect have been discussed in [114 I12j . The 
above are all examples of stationary world-lines, for which the response function is constant. 
Our main goal is to get some insight into properties of the response function for more general, 
non-stationary world- lines. 

All along we assume that the interaction of the detector with the relevant field is ever present 
and constant with time. Otherwise the response function will be affected. For example even 
inertial detectors will be characterized by a non-zero response function if the interaction strength 
is variable. 

The metric is (+,—,—,—) and, except when expressed explicitly, the value of the constants 
(c, h,G,kB) is taken as unity. 

Section 2 is dedicated to the discussion of some general properties of the response func- 
tion including our proposed regularization. This is then applied in section 3 to derive various 
numerical and analytical results for non-stationary trajectories. For example, we give integral 
representations of the transition rate as function of the ratios E/a and E/T in 1 + 1 and 2 + 1 
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dimensions. We also constrain the asymptotic behavior of the response function in the high E 
regime. Section 4 summarizes our conclusions. 

2 The response function 
2.1 Preliminaries 

We start with some well documented preliminaries. One considers a massless, chargeless scalar 
field &(t,x). In a flat Minkowski space-time one can use particle states associated with inertial 
observers and decompose the field operator into these states according to the standard definition: 

*(*,*)= 1 («* + a\e^) , (1) 

JJJ-oc ^2(27r) 3 p|| 

where the annihilation and creation operators satisfy the usual commutation rules [a%, at] = 

5(k — k'). We couple this field to a point-like moving detector constrained to follow the time-like 
world-line x^(t) = [t(r),x(r)], where r is the proper-time along the trajectory. 

Following the two-level-detector paradigm, the detector's intrinsic structure is described by 
two levels: the ground level | j) and an excited one 1 1), separated by an energy gap E in the 
proper frame. It is standard to take the interaction Hamiltonian for the model, in the interaction 
picture, as 



H mt (r)=9 (e lET \])(i \ + e"^ |i)(Tl) *W(t)] , (2) 

where g is the coupling constant. By this expression, H,i nt can be interpreted as a sort of 
boundary condition imposed on the field along the world- line x^(t). The state of the system is 
the product of the field's particle content and of the detector's state. Denoting the field vacuum 
by |0), the ground level of the entire system is |0, j). Suppose the system is in this state at time 
ro: |*(to)) = |0, 1). In the absence of interaction this is an eigenstate of the free Hamiltonian 
and so remain stable. When the interaction is on, the field vacuum becomes unstable and gets 
populated by quanta while the detector makes transitions. Starting from the system vacuum at 
tq , one gets, at a later time r, a system state 



Tdr' H int (r) 

J TO 



Mr))=Te Jr a |o,|), (3) 

where T stands for time ordering. The probability P+(r) to find the detector in the upper 

level at r (independently of the field particle content) is given by the expectation value of the 

projector n + = |T)(t| i n the above mentioned state. Another useful quantity-the subject of 

d-P 

our discussion-is the transition rate (response function) defined as R+(t) = ^ . Retaining the 
lowest (second) order contribution in the interaction strength g, we have for t > tq 

P+(T) 55 (*(T)|t><t|*(T)> (4) 

= g 2 I dn I &T 2 e- iE ^-^W{T 1 ,T 2 ) 

J TO J T() 

= 2g 2 Re dn / dr 2 e~ iEr2 W(n, n - r 2 ) (5) 
Jr Jo 

R + (t) = 2g 2 Re / dr' e~ iEr ' W{t, t - r') , (6) 
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where 



W(n,T 2 ) = (0|$^(ri)]$[^(r 2 )]|0) (7) 

is the Wightman two-point (vacuum) function evaluated along the trajectory, and use was made 
of its property: W(ti,7^)* = W(t2,t\). The Wightman function can be written explicitly in 
terms of the trajectory: 

W(t u t 2 ) = e -M^(n)-^(r 2 )] (8) 

2(2tt) JJJ-oc 



The probability P_ and the corresponding response function i?_ related to the de-excitation 
process are obtained by replacing E by — E in the expressions above. If a steady-state is reached, 
the ratio of level populations is n + /ri_ = 

Equations© and ([6]) exhibit causality explicitly: to get these quantities at time r we only 
have to know the trajectory prior to this time. 

The quantities we deal with here enjoy an obvious and useful scaling relation stemming from 
invariance under change in the units of length: If two world-lines, x^{t) and x^(f) are related 
by ^(f) = Xx fM (f/X) then their respective response functions are related by 



R+(f, E) = A _1 J R + (f/A, \E) . 



(9) 



The intrinsic, Poincare invariant, quantities that describe a world-line x^(t) in 3 + 1 space- 
time can be chosen as the acceleration a(r), the torsion T(r), and the hyper-torsion H(t), 
defined as 



a(r) 
T(t) 
H{t) 



a -1 [a 4 — a 2 



i T -2 e ^jS ± ■ 



apa^as , 



(10) 

(11) 
(12) 



where ' = a M = x^(r), and e a P~<& is the Levi-Civita anti-symmetric tensor (compare Letaw 
|13|). A world-line for which the three invariants are r independent is called stationary. For 
such world- lines all points are equivalent and, clearly, R+(t) is constant. 

The response function, being itself Poincare invariant, can be written as a functional of these 
[instead of one of the frame specific x^(t)]. The above scaling relation implies for f(r) = a, H, T 



/(f) = A-V(f/A) 



so that 



R + [r, E; <i(t),T(t),H(t)] = A" 1 R + [A -1 r, XE; Aa(Ar), AT(Ar), XH(Xt) 



(13) 



(14) 



(We use the convention whereby a dependent variable is a functional of the function variables 
appearing after the semicolon.) This can be used to reduce the number of variables. For instance 
if we take A = E^ 1 we get 



R + [r,E;a(T),T( T ),H(T)] = E R A 
where te = r/E. 



Et: 



a(r E ) T(t e ) H{t e ) 



E 



E 



E 



(15) 
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2.2 Regularization of the Wightman function 



As expressed in Eq.([8|), the Wightman function is not well defined since the integrand behaves 
as k times an oscillating term in the UV limit. Various regularization procedures can be found 
in the literature. For instance, one can add a small imaginary part —ie term to t(ri) — t(r 2 ), 
as in Birrell and DaviesfH]. This deformation multiplies the integrand in Eq.(|8|) by e~ ke , which 
ensures convergence. This gives a Wightman function of the form 

W bd (t 1 ,t 2 ) = --^ 1 (16) 

47H \t(Ti) - t(T 2 ) - ie\ Z - [x[Tl) - x(t 2 )\ Z 

This procedure has been criticized by Takagi[l| and by Schlicht[5] who showed that it leads to 
non-physical results even for some of the most simple examples. Takagi proposed a regularization 
based on endowing the detector with a small, finite extent, which is then described with the 
help of Fermi- Walker coordinates. Schlicht[5j improved Takagi's procedure, retaining causality 
explicitly in each step. This regularization is implemented by replacing 

where e is now a measure of the proper spatial extension of the detector. Since k^ is light-like 
and is future oriented and time-like, k^x^^r) > ensures the integrability for large \k\. The 
resulting expression for the Wightman function is 



WS( - Tu T2) 4vr 2 {x»(n) - x»(t 2 ) - ie[x»(n) + x^(r 2 )]} 2 ' 

This procedure has been extended to various space-time topologies [6] as well as for more general 
kinds of spatial extents [7] • Without the regulator (i.e., putting e = 0) these two formulations of 
the Wightman function are inversely proportional to the square of geodesic distance 

d 2 (rur 2 ) = ^[^(n) - x» (r 2 )][x» (n) - x v {t 2 )} . (19) 



2.2.1 Proper-time regularization 

Motivated by several requirements that we deem desirable in the Wightman function, we would 
like to propose another regularization. The Wightman function always appears as a weight 
function in integrations that result in observable quantities such as the transition rate, the flux 
emitted [I~5]. the conditional energy[l6], the flux correlation function[T7|, etc. Hence, W(t,t — t') 
should be viewed in Eq.([U]) as a distribution of the running variable t' . In this respect, the 
analytic properties of W is of key importance. More specifically, one needs to know the poles 
distribution of the Wightman function and their relative weights in order to express analytically 
quantities such as the transition rate. 

As we previously noticed, apart from the coincidence point r' = where it diverges, the 
Wightman function is inversely proportional to the square of geodesic distance. This means that 
the analytic properties of W (its poles) come directly from the zeros of the geodesic distance in 
the complex plane. We see then that the Wightman function formally possesses only one real 
pole at t 1 = and (possibly) an infinity of poles symmetrically distributed above and below 
the real axis, since the geodesic distance is real. The regularization is only needed to cure the 
divergence at r' = 0. More exactly, the regularization's role is to isolate this pole such that 
it does not contribute to -R+, only to As a consequence, one finds that the difference 

between the excitation and de-excitation response functions, for every trajectory, is precisely 
the contribution of this pole (see next section). 
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From the last remarks, one deduces that it is desirable in a regularization that it displace the 
t' = pole above the real axis such it does not contribute to and preserve the distribution 
of other poles in. The above two procedures satisfy the first requirement but not the second. A 
way to shift the zero pole while preserving the characteristics of the world-line in a consistent 
way is to introduce the regulator by replacing x^(r) with x^(t ± ie). The wave- function in the 
integrand of Eq. ([8|) becomes 

e -ik^{T) _^ e -ik^{T-it) ^ ^20) 

which gives the regularized Wightman function 

WP\n,r 2 ) = —-r-r-f . / ~, M ., ]2 . (21) 

4"7T Z [X^[Ti — It) — X^{T2 + ie)\ z 



Schlicht's Wightman function is the first-order expansion of our Eq. (|2ip in e. This means that 
the two regularizations have the same e — > limit and should give the same results for observables 
in this limit. The superiority of one choice over the other may, however, show up when we want 
to (or have to) calculate with a finite value of e, e.g. for computational convenience. The non- 
vanishing of e matters only near the coincidence point r' = where it appears in the denominator 
in Eqs.(HH|) and ([21]) respectively as -4e 2 and -4e 2 + 4e 4 a 2 (r)/3 + C(e 6 ). Therefore, Schlicht's 
Taylor expansion agrees with ours (barring the yet higher order terms) as long as e C l/|a(r)|. 

Consider, for example the canonical linear motion with constant acceleration ao, for which 
all the expressions can be calculated exactly for finite e. The motion is hyperbolic, with 



x"(r) 



sinh(a r) cosh(a r) Q Q 



a a 
Schlicht's regularization gives for finite e 



(22) 



a 2 



Wg c (r, t-t') = -— ^ [sinh (a r' /2) - m ecosh (a r' /2)] , (23) 



16tH 

whereas our proposed proper-time regularization Eq. (121j) . gives 



W? cc (t,t- t>) = S mh- 2 (G (r' - 2ie)/2) , (24) 



The corresponding transition rates are 



— arctan(ean) 

an v 1 



e a o 



= R + X - i Z ,2,2 (25) 
R p +(e) = R u + x ~ ~ (26) 

*? - 1ST (27) 

e a o — 1 

where BH is the Unruh [3] expression, which is the e = limit of the two regularization schemes 
([x] is the integer part of x). To obtain our result, we have used the pole decomposition 
l/sinh 2 (x) = Xm=-ooV (x + inir) 2 and integrated by the method of residues. As expected, 
expressions (|25p and ([26j) agree to first order in a^e for which they both differ from the Unruh 
expression by a factor e 2Ee . So, both schemes still require eE -C 1. However, for large values of 
eao Schlicht's formulation leads to a vanishing transition rate whereas ours always differs from 
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the correct expression by a finite factor between 1 and e 27TE / a °, which is near one for E -C uq 
for arbitrary values of e. 

For world-lines with bounded values of the acceleration we can choose a value of e that is much 
smaller than l/o at all times, and then the two regularizations should give very similar results. 
The situation is, however, different for world-lines with unbounded acceleration. Consider as an 
example the tractable case of a linear space trajectory with a(r) = — a/r, — oo < r < 0. For 
a/1 this is gotten, for example, for the world-line 



\l+a 



\l-a 



for a 



l + a 1 
1 a representative world-line is 



a 



H) 



l+a 



H) 



l-a 



a 



a 



,0,0 



(28) 



s"(t) = 

(i?+(r) does not depend on k, as different values of k are related by a Poincare transforma- 
tion. This trajectory is particularly interesting because when followed by a moving mirror it 
generates a constant thermal flux oc k 2 with temperature oc n [18J.) 

For these world-lines, the scaling law, for a finite value of e, becomes (in both regular- 
ization schemes) R + (e/X, r/A, a, EX)/X = R+(e,r,a, E). Choosing A = 1/a we can write 
R+(e,T,a,E) = aR+(ae,a,E/a). The question we ask at this junction is: "how well does 
the expression for finite e approximate the correct (e = 0) value?". We present in Fig. ([I]) the 
departure factor R+(e, t, a, E)/R+(Q, r, a, E) = R + (ae,a,E/a)/R + (0,a,E/a) as a function of 
x = ae, for several values of the pair a, E/a using Schlicht's regularization, as well as ours. We 
see that, in line with what we saw in the constant acceleration case, at least for low energies the 
proper-time regularization works well for values of ea even an order larger than for Schlicht's 
regularization. 



log(- 



-TK I 



2k 



T K log( — TK,) 



T K 



2k 



-,0,0 



(29) 




Figure 1: The ratio R+(e,r, a, E)/ R+(0,T,a, E) — R+(ae, a, E /a)/ R+(0, a, E/a) as a function of ae 
for 2nE/a — 10 -2 (left panel) and 2nE/a — 1 (right panel) and for two values of a = 5, 100, as marked. 
Continuous lines correspond to the proper-time regularization, dashed curves to Schlicht's regularization. 



2.2.2 The subtracted Wightman function 

The troublesome singularity in the Wightman function that we wish to regularize appears in 
the same way (same pole position and residuum) in all trajectories including the inertial ones 
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for which the transition rate vanishes. Being time-like, all trajectories have for small r lapses: 
d 2 (r,T — t 1 ) = t' 2 + C(r /4 ). (The r' 3 term vanishes identically for all world-lines.) So there is 
some / such that 

d 2 (T,T-T')=T' 2 +T' 4 f(T,T-T'), (30) 

where / is regular at r' = ; in fact /(r, r) = a(r) 2 /12. We can then write using Eq. (|2ip 

S* r A , e-' Er ' a 2 r. , cos(Et')/(t,t-t') 
R + (r) = -— 2 Re^ dr + — 2 ^ ir ^-^——^ , (31) 

where we put e = in the second term, which is now without the pole at r' = 0. The first 
integral vanishes for positive E and e, as we close the integration contour around the lower 
complex plane. (In fact, this term is just the response function for an inertial world-line.) For 
the de-excitation rate we have 

R-ir) = 9 ^ + R + (r). (32) 

The same pole free expression was used in |10t [T3] for stationary orbits, and was generalized to 
non-stationary world-lines, though with some restrictions, in [7j, using Schlicht's regularization. 
These authors extended the applicability of the pole-free expression to finite time interaction 
[Eq.(6.1) in [7]]. We note that this finite interaction time expression follows straightforwardly 
in the proper-time regularization scheme [Eqs . (i30|) (f3TI) ] applied to the general response function 
Eq.©. It turns out that the pole- free expression is not necessarily more useful, in numerical 
determinations of the response function, than the finite e regularizations. The reason is that the 
function / in the integrand is defined as the difference between two diverging expressions and 
its numerical evaluation requires, in the end, some sort of regularization. 



3 The transition rate— some derived properties 

In this section we discuss two issues regarding the response function of general world- lines. The 
first is its behavior for high detector energy gaps for stationary and non-stationary trajectories. 
The second issue concerns the applicability of the "quasi-stationary" approximation, whereby 
the response function is represented by that of a stationary world-line with the instantaneous 
values of acceleration, torsion, and hyper-torsion. 

3.1 High energy behavior 

To study the high energy behavior of the response function we use the subtracted Wightman 
function Eq. ([3Tj) . One can then write: 



R-\ 

where 



/•CO 

/ dr cos(Et')F t (t'), (33) 
Jo 



FAt) ~ 2^l + T' 2 f(r,r-T'Y (34) 
Following an idea in [7] we integrate by parts 2N times to get for any N 
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*f(r) = £ (^) n ^ 2n " 1} (0) + (^) 7V fdr' co S (Et>)F? n Ht'), (35) 

where we used the fact that Fj n ^(r' — ► oo) = since F t {t') decreases at least as rapidly as t /_2 . 
Using the Riemann-Lebesgue lemma [T9], the second term is a 0(x~ 2N ^. So, Eq.([35]) can be 
used to pin down the asymptotic behavior of the response function for E —* 00. This we proceed 
to do below, first for stationary trajectories, then for non-stationary ones. 

The first three odd derivatives of F T are useful and are given in terms of the world-line 
derivatives as: 



£)"V(0> = (30) 



9 



_ M 3) (0) = 1 1 1 , (37) 

2vr 2 / 12 20 30 60 60 v ; 

(5) a^a a 3 aT 2 13ad 3 a 4 rT 

2^2y r W ~~2T 9 63 HT 

5a 2 TT 5aaf 2 53a 2 da 5ddT 2 
+ 42 + 42 + 126 + 84 

5adTT 5adrr 5a 2 ff aV 3 ) 
+ 42 + 42 + 84 + 14 

5aa( 3 )T 2 5da( 3 ) 5a 2 Tr( 3 ) 5dq( 4 ) 

+ 252 42" + 252 72" 

(5) ^(4) (5) 

+ — ^— . (38) 



56 168 



3.1.1 Stationary trajectories 



For stationary world-lines the geodesic distance, Eq. (|19p . is a function of the lapse r' = t\ — T2; 
so, by Eq.©, the transition rate is constant (provided the interaction strength is constant): 
We can write with the proper-time regularization, using W(t, t — r') =>- W stat (T r — ie), and the 
general property W stat (T')* = W stat {-r') 

/+00 
dr' e~ iET ' W stat (r' -ie) , (39) 
-00 

These world-lines were classified and analyzed for example in [T3] . They can be grouped into 
six classes, equivalent to the orbits of time-like Killing vector fields. As said above, they are all 
characterized by constant values of a, T, and H. For a = T = H = we get inertial world-lines; 
linear, uniformly accelerated ones have T = H = 0, a / 0; and circular, constant speed ones 
have (a < T, H = 0) (H = for planar trajectories). 

The general property d 2 (ri,T2) = d 2 (r2,ri) leads, for stationary trajectories, to F t (—t') = 
F t (t'). Thus, all the odd derivatives of F T at r 1 = vanish, and with them the first sum in 
Eq. (|35p . for any N. We conclude then that for all stationary trajectories i?+ vanishes faster 
than any power of E in the high energy limit. This is clearly in line with the known thermality 
of the response function for uniformly accelerated trajectories (the original Unruh case Eq. (|27|) ). 
We discuss other examples below. 
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Note first that from the general scaling law (|15p . we can reduce by one the number of 
independent variables in the stationary case. For example we have in this case: 



R+(E,a,T,H) = E R + (a/E, T/E, H/E) . 



(40) 



For general planar trajectories (up to a Poincare transformation) the hyper-torsion vanishes, 
and the stationary trajectories are characterized by the acceleration and the torsion. Circular 
trajectories, which are defined by T > a, can also be defined by some orbital radius R and a 
frequency Q in the frame in which the center of the orbit is at rest, such that a = ^y 2 Q 2 R and 
T = ja/y/ 7 2 — 1, where 7 = (1 — R 2 $l 2 )~ l l 2 . T < a defines hyperbolically unbound trajectories, 
and T = a, which is the the 7 — ► 00 limit of the circular case, corresponds to the so called "cusp 
motion" . It is realized, for example, by the world-line 



J cusp 



T + 



a 2 r 3 



ar 2 a 2 r 3 



6 



-,o 



A circular world-line is realized, for example, by: 



X circ( T ) 



Tt a 

, — ^ COS (UJT) 

U) L0 Z 



a 

or 



sin (lot) , 



(41) 



(42) 



where u = (T 2 — a 2 ) 1 / 2 (for circular orbits u = 7^). The geodesic distances for these three 
types of world-lines is given by 

4a 2 



d 2 (r,r 



T 2 



2 12 

LO T 



T 2 



sin 



K/2 



(43) 



and is applicable also for imaginary (hyperbolically unbound case), and vanishing (cusp case) 

LO. 

The response function for the three stationary cases is of the form E times a function of 
v = a/E and fi = T/E. Defining y = Er/2 and k = lo/E = (n 2 — v 2 ) 1 / 2 we obtain with the 
proper-time regularization: 

9 2 E 



Ra 



4vr 2 



dy 



_ K 4 e -2iy 



K 2 /J, 2 (y — it) 2 — v 2 sin 2 [«;(y — it) 



An analytic expression is known for the cusp motion [13]: 



R cusp 



g 2 a e - 2 ^ E ' a 



8vr\/3 



(44) 



(45) 



which, again, decreases with E faster than any power, as we found for all stationary trajectories. 
This is also the response function for the circular case in the limit of infinite 7, which was derived 
before by Bell and Leinaas[8] and by Unruh[9] . (They proceeded by reducing the circular geodesic 
distance Eq. (|43p to the contribution of its first three zeros, integrating the Wightman function 
by the method of residues. It can easily be understood that this gives the same result as using 
the cusp geodesic distance Eq. (j4"3|) for lo = 0.) This infinite 7 "circular Unruh effect" is quasi 
thermal in the sense that it can be described by a "temperature" that varies slowly with E, 
monotonically from T < = at low energies to T > = at high energies. Letaw and 
Pfautsch[10j used the subtracted Wightman function to numerically calculate the behavior of 
the transition rate for finite a and T. It is of interest to have an analytic expression for R" rc 
for finite values of 7. We derive the first order correction in 7 -2 . For the geodesic distances one 
finds 



df. 



d 2 \l 

""cusp L x 



4 4 

a t 



360 + 30a 2 r 2 7 2 



+ u(ar)0(r% 



(46) 
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where d 2 usp is the expression for the cusp world- line [Eq. ([15]) with u = 0] ; From this one deduces 
for fixed E/a in the limit of large 7 



i?^(£,a, 7 )=iC P (£,a) 



1 + ^(^- 1 ) + ^ /0)0(r " ) 



(47) 



We verified this with numerical results for i? + . Using this and the expression for i?_ we see that 

"7 



the effective temperature in the limit E — > is modified to T< ~ T < (1 — I/57 2 ). We do not 
have an analytic expression for fixed 7 and high E. 



3.1.2 Non-stationary trajectories 

We start again with the subtracted Wightman function as in Eq. (135p . We noted that for a 
stationary trajectory all the odd derivatives vanish because the expression for the geodesic 
distance from r to t — t' is symmetric in r' at r' = 0. The opposite is also true: if there is such a 
symmetry for all values of r the trajectory is stationary. Thus, if the trajectory is not stationary 
there will be r values (generically, all values of r) where at least some of the odd derivatives 
do not vanish. For such r the response function decays only as a power-law at high energies, 
and thus, in particular, the spectrum is not thermal. For example, if d 7^ the dominant high 
energy behavior, from Eq. (l36p . is 

R + (r) « <(r) = ^9- (48) 

If the acceleration is constant on a world- line but the torsion is not, the dominant contribution 
is, generically, proportional to the derivative of the torsion T(r). From Eq. (|37p we get 

, . g 2 a?TT , . 

(Note that if a vanishes at an isolated point, the first term in Eq. (|37p gives a finite contribution 
to the response function at that time.) If T is also constant the next order will dominate and 
it contains the derivative of the hyper-torsion through x^^xJ^ oc HH. If this too vanishes 
the trajectory is a stationary one. Conversely, if the world-line is non-stationary, so either a, T, 
or H do not vanish identically, at least one of the first three orders doesn't vanish identically. 
We conclude that generically the high-energy behavior of the transition rate for non-stationary 
trajectories is given by R+{t) ~ g 2 a(r)E~ 2n , where a(r) contains the acceleration, the torsion 
or the hyper-torsion and their derivatives and where n = 1, 2 or 3. 

We checked numerically the validity of the first order approximation for the linear tra- 
jectory Eqs.(|28j)(|29]), having a(r) = -a/r. In Fig.© we plot the ratio £ = R+/R+, cal- 
culated numerically with the proper-time regularization. Using the scaling law we can write 
R+(r,a, E) = ER + (E / a,a). Since R+ is obtained from a Taylor expansion of R + , it also sat- 
isfies this scaling. Thus, the ratio £ = i?+/i? p is a function of a and E/a. Fig. ([2]) shows this 
ratio as a function of E/a for different values of a. We see that indeed Eq. (|48p is a very good 
approximation for the response function at large values of E/a. For a larger than about 5 the 
power-law behavior starts at the same value of E/a regardless of the value of a. In the large-a- 
low-E/a regime £ can, in fact, be approximated analytically very well using the quasi-stationary 
approximation-see section 3.2 and Fig. ([5]), so the behavior displayed in Fig. ([2]) for this regime 
can be reproduced analytically. 

In the above example, once E/a 3> 1, the dominant power in E/a is a good description of 
the response function for all values of a > 1. This can be traced to the fact that in this case 
E~ 2 fP(0)/F!t(0) = 0[(a/E) 2 } (T = for the present, linear case). In general, the contribution 
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Figure 2: The ratio £ = R+/R^_ for the linear trajectory Eq. (|28[l as a function of E/a, for different 
values of a (as marked). 
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Figure 3: The response function for the linear trajectory with acceleration a(r) = 2 + tanhr for E = 20 
(g — 1). The result is given in units of 10 -5 . The dashed vertical line shows the maximum of BF, 
to = In [(2 + y/7)/3] /2. 



of the higher orders hinges also on other characteristics of the world-line as reflected, e.g., in the 
additional terms in F T (0). To better demonstrate this we considered another example: Take 
a linear trajectory with acceleration law a(r) = 2 + tanhr in some arbitrary units (r in the 
inverse of the same unit); a is bound between 1 in the far past and 3 in the far future. Such an 
acceleration law is gotten for example for the linear world-line: 



+ — + — - + arctan(e~ T ), e~ T + — + — — arctan(e~ r ), 0, 



(50) 



At large negative times the detector has a nearly constant acceleration with the characteristic 
time for acceleration change t± r = a/a ~ (3/4)e~ 2r 3> t a = 1/a ~ 1. The first term in the 
expression for (0) dominates again and then we expect that for S>1 the first order power- 
law will be a good description. As r approaches zero the additional terms in the expression 
for Ft (0) become important. We first show in Fig. ([3]) the computed transition rate itself 
for a single value of the energy E = 20 (in the same units). It follows tx a(T)a(T) = 
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Figure 4: The ratio £ = R+/R^_ for the linear trajectory Eq. (|50[) with acceleration a(r) = 2 + tanhr 
and for different values of the energy. 



(2 + tanhr)/ cosh 2 r rather closely (as expected since E S> a(r) for all r). The latter is maximal 
at t = In [(2 + \/7)/3] /2 at which it takes the value 2(10 + 7^/7)/27 ps 2.11. 

In Fig.Q we plot the ratio £ = R+/R+ for different values of £7. For i£ larger than about 
5 all the curves cross a value of 1 at about the same value of r ~ —0.906313, which is a zero of 
Fj 3) (0). 

This is a reminder that all our results are valid only to the lowest order in g 2 . So, for 
example, if we speak of asymptotic behavior it is to be understood in the following sense: Given 
a high enough energy, E, we assume that g is small enough for the lowest order to dominate 
at E. We have no reason to think that our results on the asymptotic behavior are invalidated 
by higher order corrections. However, even if they are, our claim is even strengthened. Take 
the non-stationary case. If higher order terms decay even faster with E, then obviously they 
do not affect our results. If, in some case, higher order terms decay even slower than the first 
order behavior we found, then this even strengthen our conclusion regarding the non-thermal 
behavior and the power-law decay. In the stationary case where we found the first order response 
function to decay faster than any power, it might in principle happen that higher order terms 
decay slower. However, this is known not to be the case at least for the Unruh case where R + 
is known to all orders. 

3.2 The quasi-stationary approximation 

One might think that if the world-line is characterized by a, T, and H that vary slowly enough, 
the response function is well approximated by the expression 

R+[T,E;a(T),T(T),H(r)] ^ R q + [E,a(r),T(r),H(r)} = E R + 

Here, the value of i?+ at r, which is a functional of a(r), T(t),H(t), is approximated by 
its value for a stationary world-line with the instantaneous values of these parameters. For 
example, Letaw |13| remarks on linear trajectories: "if a detector was moving with a very slowly 
increasing linear acceleration one would expect the spectrum to be Planckian, though with a 
slowly increasing temperature". Our analysis of the high energy behavior of R + shows that 
this cannot be correct across the spectrum. For any given non-stationary world-line, no matter 
how slow the variations in a, H, T are, the quasi-stationary approximation has to break down 
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E 1 E 



E 



(51) 
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at high enough energies as the spectrum becomes power-law as opposed to that of stationary 
world-lines. However, at low enough energies, slowness of variation of the three invariants might 
suffice to validate the quasi-stationary approximation. It is difficult to give general criteria for 
this so we treat numerically several linear motions. In this case we write 

R + (r,E) = R%{t,E) + D{t,E) , (52) 

R%(r,E) = ^ — 4 , (53) 

e K T >I — 1 

where R\_ is the stationary (Unruh) expression for the response function. We wish to know 
under what conditions \D{t,E)\ <C R q + {T, E). As we said above, one necessary condition is that 
E is small enough, and it is convenient to express this smallness in terms of E/a (provided 
locally). Another condition it that the acceleration varies slowly enough on the world- line: the 
time scale of variation of a: tf r = a/a, has to be long compared with t a = 1/a. In other words 
a /a? <C 1. In general, tt r is not the only measure of the orbital time variations, but in the 
examples that we consider below it is a good single measure. We considered, numerically again, 
the two families of world- lines discussed before. In the first instance we looked at the world-line 
with a(r) = —a/r. In this case, the time ratio a 2 /a = a is constant along the world- line, and 
it also represents, up to some numerical factor the time ratio for all the possible time scales 
defined by the trajectory. We show in Fig.(|5]) the ratio £ = R+fR 9 ^ as a function of E/a for 
several values of the parameter a. Because of the scaling properties mentioned above, C is a 
function of only these two variables, and is given by: 

e 2nE/a _ ^ i-oo ^ f a 2 - I 1 ) 

C = ^r; / d-r' cos (ar E/a) I ; — ; -— — - -\ -r > . (54) 

S iraE/a J V 1 ' \ [1 - (1 + r') 1 ^] [1 - (1 + r') 1+a ] r' 2 J y 1 
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Figure 5: The ratio £ = R+/ ' R q + as a function of E/a for different values of a (as marked) for the linear 
trajectory with acceleration a(r) = —a/r. 



We see that indeed the quasi-stationary approximation is good (£ ~ 1) when E/a is small 
and the orbital variations are slow (high a), but it breaks down for all values of a plotted for 
large E/a, and for all values of E if a is small. The parameter £ has a well defined low energy 
limit, which for high a behaves, according to the numerical results, as C(E = 0,a) ~ 1 — a . 

Next, we reconsidered the linear trajectory Eq. (|5U|) with acceleration a(r) = 2+tanhr. Since 
the acceleration is bounded between 1 and 3, the condition E ^ a is tantamount to E <g 1. The 
ratio of time scales a/a 2 is small for large, positive and negative, values of r where it decreases 
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Figure 6: The ratio C, = i?+ /-Ri for different values of the energy for the linear trajectory with a(r) 
2 + tanhr. 



exponentially, but is of order one near r = 0. In Fig. ([6]) we plot £ as a function of r for different 
values of E. For E = 2-which corresponds to E/a ~ 2 for large negative r, and -E/a ~ 2/3 for 
high positive r-we see that C departs from one appreciably even when the acceleration varies 
very slowly. Also for very small values of E the quasi-stationary approximation is valid only 
when the acceleration varies on time scales longer than 1/a (i.e., away from r = 0). 



Finally, we consider a linear trajectory with an acceleration law a(r) 
0). It is realized, e.g., by the world-line 



-2/3 



(_ T )-l/3 (r < 



z"(r) 



TO 



dx sinhx 2 - ^sinh^ 2 ^ , ~to^ (J dx cosh x 2 - A cosh A 2 ^j , 0, 



,(55) 



1/3 



where ^4 = y | ^-^J . This acceleration law might be of interest because it reproduces the 

time dependence of the surface acceleration of an evaporating black hole (in which case r is the 
proper-time of a static observer at infinity). Accelerating detectors and black hole are related 
(e.g., An analogy can be drawn between the Unruh thermal bath for a uniform acceleration 
a, and Hawking radiation from a constant mass black hole with surface gravity k = 1/4M, as 
the temperatures for the two cases are given respectively by 



ha hn 
kbTu = - — , k B T H = - — 
Zttc Zirc 



(56) 



If we apply the semiclassical treatment down to complete evaporation through a massless scalar 
field, and take the time r = at this event, one finds that k = T p 2 ^ 3 (— r) -1 / 3 , where Tp = 
tpi/\/80ir and tpi is the Planck time (the numerical factor relating Tp and tpi varies when 
considering other fields according to their mass, spin, etc). 

In discussing black hole evaporation one usually assumes the analogue of our quasi-stationary 
approximation; i.e., that at any given time a black hole radiates as a black body corresponding 
to its momentary mass. We set to check to what extent such an approximation is justified in our 
mock-evaporation cas^ll Of course, our chosen acceleration law is based on the quasi-stationary 



x The process of Hawking radiation has been proven to keep its thermal character until the mass of the black-hole 
reaches O(Mpi) (see [2U| for a heuristic picture, [3T| for a direct numerical treatment and [22] for a semiclassical 
back-reaction analysis). Our analogy is an alternative, heuristic check. 
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approximation in deducing the time dependence of the mass. It is also clear that the analogy 
becomes meaningless beyond about r = —Tp. But we still consider the analogy at all times 
as it may be of interest for its own sake: in general our parameter tq need not be related to 
the Planck time. Again, we consider the parameter £ = R+/R+ as a measure of the departure 
from quasi-stationarity. The scaling law for R + tells us that this quantity can be written as a 
function of two parameters; for example £ = ((t/tq, E/a), or alternatively ( = C(t/tq,Etq). 
The first helps to follow the time dependence of R+ at a fixed value of E/a, which tells us where 
on the black body spectrum we are; the second representation follows the time evolution for a 
fixed value of the energy. The time variable —t/tq = (3d/a 2 ) _3//2 is also the parameter we use 
to measures the fastness of the variation of a. 
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Figure 7: The ratio £ = Rj r /R q + as a function of t/tq for different values of E/a (as marked) for the 
linear trajectory with acceleration a(r) — r 2 ^ 3 (— r) -1 ^ 3 . 




We plot the above two representations of £ (computed with the proper-time regularization) 
in Figs. ([7]) and © respectively. In the first ( is plotted vs. t/tq for different values of E/a and 
in the second for different values of Etq. (For values of these parameters below the minimum 
ones plotted the curves remain the same.) We see that, in line with our previous results, 
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for high enough energies, or rapid enough variations in the acceleration, the quasi-stationary 
approximation breaks down appreciably. For —t/tq > 100, where a /a 2 is about 0.015, the 
quasi-stationary approximation is valid for all energies around and below the black body peak 
at E/a = 1/2-7T. In the black hole analogy this value of t/tq would correspond to a black 
hole mass of ~ 0.07Mp; (Mpi is the Planck mass). So our results would indicate that the 
quasi-stationary approximation for black holes holds in the whole classical regime. 

One can show that, in the limit r — » 0, -R+ tends to an energy-dependent finite value RP^(E) 
(unlike the a(r) = — a/r case, where i?+ diverges in this limit). We found numerically that 

HP + ( E ) (2/157t)( 3 2 /t ). 

We see then that at least for the above three examples a /a 2 <C 1 is indeed a sufficient 
condition for the validity of the quasi-stationary approximation, provided E/a < 1. 

4 Conclusion 

We considered several issues concerning the behavior of the transition rate for a moving two- 
level detector in the vacuum. We introduced the Feynmann prescription for regularizing the 
Wightman function in order to secure causality for the response function of a two-level detector. 
We then use this to derive certain properties of the response function. We showed that for 
stationary world-lines the response function decreases faster than any power-law in the high 
energy limit, in concordance with the known exponential decrease for a couple of known special 
cases. In contrast, non-stationary motions are characterized by a power-law decline of the 
response function at high detector energies: Generically, the response function at any given 
time vanishes as E~ 2 at high energies; in special cases the decline is faster: either as E~ 4 or as 
-E -6 . We also considered some aspects of the applicability of the quasi-stationary approximation, 
whereby the response function is approximated by it value for a stationary trajectory with the 
instantaneous values of the invariant parameters a, T, H. Our result on the high energy 
behavior precludes this approximation for high energies (generically 1 < E/a), but in the low 
energy regime we find for several examples studied numerically that the approximation is good 
when a/a 2 <C 1. 
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